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This paper is presented in two distinct parts. 

Part 1 begins by giving a number of examples of random surfaces and summarizing the 
various applications of random surface theory,, including applications to 

(a) electromagnetic scattering from such surfaces, 

(b) the information content of photographs, maps, etc. interpreted as random surfaces, 
and implications for bandwidth compression. 

This part of the paper will also summarize the major results which have previously been 
obtained, and will outline the approaches used to obtain these results. Finally, some inter- 
esting additional problems in this field are stated. 

Part 2 of the paper is devoted to new results on the properties of random spherical 
surfaces, and the application to electromagnetic scattering from such surfaces. These results 
primarily concern the statistics of the number of specular points per unit surface area. 

The detailed elaboration of the mathematical approaches to obtaining the major results 
of random surface theory, and even the exhibiting of the major formulas, would be very com- 
plex and space consuming. Therefore, insofar as previously obtained results are concerned, 
we will merely give outlines and examples, with references to the appropriate literature for 
the details. Such detailed mat hematics as appears in this paper is reserved for the new 
results presented in part 2. 

For convenience, the two parts have separate numbering for equations and references. 

Part 1 — Outline of Random Surface Theory and 

Applications 

1. Introduction 
1.1. Definition and Fxamples 

A random surface is a random process {/(/,//)} depending- (usually) on two real parameters 
./' and ?/. Although we have used the notation (x,y), these 1 do not necessarily represent 
rectilinear Cartesian coordinates. 

For example, one of the parameters may represent time: y=t. In this case, we ma} 7 be 
dealing with a random function of one spatial parameter winch is also a random function of 
time. However, the mathematical techniques are those applicable to random surfaces. 

Also, the parameters, even if they represent spatial coordinates, may represent coordinates 
other than Cartesian. For example, we may have x=d, y=<f>, where 0, <j> are spherical angular 
coordinates. 

Occasionally, we may be dealing with a random function of two spatial coordinates which 
is also a random function of time, in which case the random surface in question must be repre- 
sented as a stochastic process {f(x,y,t)} depending on three real parameters. 

Examples of random surfaces are: 

(a) The surface of the ocean or other body of water (in which / represents level with respect 
to a mean level), 

(b) a topographic map (/represents altitude), 

(c) an atmospheric pressure chart (/represents atmospheric pressure), 

(d) a photograph or television picture (/represents brightness), 



1 A portion of the work reported on in this paper was supported bv the National Aeronautics and Space Administration under Contract 
NASA 490. 
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(e) an antenna pattern in two spherical angular coordinates, and 

(f) a rough spherical surface, such as that of the moon or other planet considered as a 
whole. 

All of the above, but in particular, (a), (c), and (d), may depend not only on two spatial 
coordinates but on time as well. 

As previously stated, random functions of two parameters, one of which represents a 
spatial coordinate and the other time (as in a moving waveform) can also be considered to be 
random surfaces from the point of view of statistical treatment. 

1.2. Applications 

The most extensive treatment of random surfaces remains that of the series of papers of 
Longuet-Higgins [1952; 1956; 1957 a, b, c]. Much of this was motivated by applications to the 
surface of the sea. 

In particular, many of the properties studied are of interest for studying the electromag- 
netic (radar or optical) scattering properties of the sea. For example, the statistics of specular 
points and their motion are properties affecting the intensity and spectrum of scattered energy 
in a given direction. 

Clearly, this application to electromagnetic scattering properties is not restricted to the 
surface of the sea, but may be relevant to scattering from any object which can reasonably be 
represented as a random surface. Part 2 of this paper deals with the statistics of specular 
points on a random spherical surface, and was motivated by problems of reflection of radar 
waves from the moon. 

Another major application of random surface theory is to the information content of 
various examples of random surfaces (in particular, photographs and television pictures but 
also including, possibly, topographic maps, pressure charts, antenna patterns, etc.). 

By information content is meant how many bits of information are required to describe 
such a surface to some given degree of accuracy. This affects the storage requirement for 
storing large amounts of data of this type, as well as the bandwidth requirements for transmit- 
ting such data through a communication channel. The work of Swerling [1962] on contours of 
random surfaces was motivated by these considerations. 

If the statistics of the random surface are completely specified, the information content 
can be determined according to the standard formulas of information theory [Swerling, 1962]. 
However, bounds on the information content may be more easily obtainable by studying 
certain statistical properties of contours, such as expected length per unit area, expected number 
of separate closed contour pieces, and expected "roughness" of contours. 

The above discussion is by no means exhaustive and is in fact only a sketchy outline of 
the possible applications of random surface theory and of the particular properties of random 
surfaces which are of interest in these applications. For example, extensive study has been 
made of the statistics of the number and velocity of maxima and minima of moving waveforms, 
the velocity of the crossings of any given level, the annihilation and creation of specular points, 
etc. [Longuet-Higgins, 1952; 1956; 1957 a, b, c]. 

2. Analytical Techniques for Random Surfaces 

The statistical properties of random surfaces which have been most extensively studied 
involve surface properties which can be described in terms of the differential geometry of the 
surface — that is, in terms of the values taken on by the function /and its partial derivatives of 
various orders. 

These properties involve the probability that the function/ and its partial derivatives take 
on certain values or that the line integrals, along contours, take on certain functions of the 
partial derivatives. 

For example, an extreme point of a random surface is characterized by f x =0, j y =0. A 
contour is characterized by J=z Q . In a moving waveform j (x, t), the velocity of a point at a 
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constant level /= z x is the value of f t /f x subject to/= z x . Specular points can be defined as points 
at which the gradient has a definite value. Such examples can be multiplied indefinitely. 

Thus, the approach to deriving the desired statistics of these surface properties consists 
first in describing the phenomenon in question in terms of the values of/ and its partial deriva- 
tives, and then taking the appropriate statistical average. This will be illustrated by several 
specific examples below. 

The results depend on the joint probability density function (pdf) of / and its partial 
derivatives. It is worth noting that most of the results which have been obtained involve the 
joint pdf of the function and its partials evaluated at the same point. However, as will be 
pointed out, many interesting statistical properties involve the joint pdf of the function and its 
partials evaluated at two or more different points; these properties have not been extensively 
treated. 

As a first example, consider the statistical properties of the length of contours in a planar 
surface. Following Swerling [1962], let the intersection of the surface z=f(x, y) with a plane 
z=z be called a £ -contour. Let V/ denote the gradient, that is the vector having components 

Define 

F d (u,v)=2ii \u-v\<2 

=0it\u-v\>^ (1) 

Then, the total length L R (s ) of the z -contour in a region A* of the {x, y)-plane is 

Z*(z )=lim f \Vf\F a (J, Zo)dxdy. (2) 

d-^OJ r 

Equation (2) gives an expression for L R (z ) apart from any statistical considerations. It 
is worth noting that this expression depends on the fact that we are dealing with planar sur- 
faces rather than, say, spherical surfaces in this example. From this expression one can derive 
the desired statistical properties of L R (z Q ) in terms of those of/. 

For example, denoting expected values by E{ ), we have 

E[L R (z )]= J I -y/ul+ulWizo, u h u 2 \x } y)du l du 2 dxdy ) (3) 

where W(u 0j u u u 2 \x, y) is the joint pdf of/, f x , and/,, all evaluated at (x, y). In (3), the inner 
integral is extended over — °°<^i, u 2 <i c °. If {f(x,y)} is stationary, W is independent of 
x, y; the expected length per unit area is obtained by dropping the integral over R. 

A second example involves events taking place at isolated points rather than along a whole 
contour, namely, the expected number of horizontal tangents to a 2 -contour. Such a tangent 
occurs when/=2 > f x =0. 

The number of points within R where f=z ,f x =0 is 

T R (z ) =lim f \f x f xy -f v f xx \f d (/„ 0) F d (/, 2o ) dxdy, (4) 

d->0 J R 

and the expected number in R is 

E[Tr(z )]= I J \u 2 u z \W(z , 0, u 2 , u 3 \x, y)du 2 du z dxdy, (5) 

where W(u , u u u 2 , u 3 \x, y) is the joint pdf of f,f x ,f yf and/^ at the point x, y. 
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A third example involves line integrals of the curvature of a contour. Let a 2 -contour 
be expressed in parametric form x(s), y(s). The line integral along the contour of (0'(s)) 2 is 
a measure of the roughness of the contour, where 0(s) is the angle of a tangent to a contour. 

Then [Swerling, 1962], 

B QVWP & ] = £ j (u 2 1 +ul)^ 2 (u 2 1 u 5 -2u 1 u 2 u 4 +uiu 3 ) 2 

XW(z , u u u 2 , u 3 , ut, u 5 \x, y)du l . . . du^dxdy, (6) 
where 

W(u , u 1} u 2 , u, 3 it,*, u 5 \x,y) is the joint pdf of/, f x , f y , f xx , f xy) and j yy at x,y. 

The integral on the left side of (6) is the line integral along the ^-contour within R. 
Expression (6) arises from the fact that d'(s) can be expressed in terms of the partials of 
/ as follows: 

f\t / \ J Jxjyy ^Jxjyjxy ~t~ . 1 yjxx I /^\ 

W ~ I lv/1 3 /" {) 

These examples are sufficient to illustrate the procedure: first express the surface property 
in terms of the appropriate conditions on the function/ and its partials; then, find a nonsta- 
tistical expression for the number of points at which the property holds, the length of the 
contour, or whatever may be the precise quantity for which one wishes to derive statistics; 
this expression can then be used to derive the statistical properties of this quantity, such as 
its expected value. 

For gaussian random surfaces, the expressions in terms of the joint pdfs, W can often be 
explicitly evaluated. The probability distributions for gaussian surfaces can be specified by 
specifying the function 

<fti,v\x, y)=E[f(x+£, y+r,)J(x,y)]. (8) 

(If the process is stationary, the dependence on x, y disappears.) 

Further, all of the joint pdfs, W can be explicitly written down in terms of the function 
<p and its partial derivations with respect to £ and rj, evaluated at %=0, tj = 0. 

The reader is referred to the literature (references, Part 1) for the resulting expressions. 

3. Additional Problems 

A number of additional problems of interest remain to be treated in the field of random 
surfaces. 

One such problem is that of extending the results to nonplanar random surfaces, such as 
spherical surfaces. Part 2 of this paper gives new results on this precise problem. 

As previously noted, the specific form taken by the expressions depends on the differential 
geometry of the surface in question. 

A second problem involves deriving the higher order moments rather than just expected 
values. This problem can easily be attacked in principle by utilizing the exact (nonstatistical) 
expressions for the quantities in question e-g., (2) and (4). 

As an illustration, suppose it is desired to derive the second moment of L R (z ), the length 
of a ^-contour in R. From (2), 

[L B (z Q )Y= lim f f \Vf(x,y)\\Vj(x^y')\F d U(x,y) y 2 ]FAf(x^y') } z ]dxdydx'dy' (9) 

d, d'^>Q J R J R 

Then, E[L B (z )] 2 can be obtained directly from (9). The resulting expression involves the 
joint probability function of j{x, y),f x (x, y) , J y (x, y) , f(x f , y'),M x '> y'),fy(%', y'), for all pairs 
of points (x, y) and (x',y'). It is very difficult to evaluate explicitly even for stationary 
gaussian surfaces. 
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Similarly, higher moments may be derived by expressing the powers of L R (z ) as multiple 
integrals and taking expected values. 

The same remarks apply to other quantities for which the higher moments may be desired. 

A third category of problems of significance is the further study of the contours regarded 
as random processes. A connected piece of a contour can be expressed in the parametric form 
x(s), y(s) where s = arc length. The derivatives x'(s), y'(s) are related; in fact, y'(s) is 
determined up to sign by x'(s). The statistical properties of x'(s) or y'(s) regarded as random 
processes are of interest. Alternatively, one could study 6 (s), 6= angle of tangent to a contour, 
regarded as a random process; 6 (s) is a function of x'(s) and y'(s). 

Swerling [1962] derives an expression for the expected value of the line integral along a 
contour of [0'(.s*)] 2 - Similarly, expressions can be derived for the expected value of the line 
integral of [x"(s)] 2 or [y"(s)] 2 . It is pointed out by Swerling [1962] that these expected values 
can be used to define quantities analogous to the second moment of the spectral density function 
of 0(s) or z'(s), assuming that such spectral density functions could be defined rigorously. 

In exactly the same way, expected values of line integrals of the squares of higher derivatives 
of 6 (s) or x'(s) can be obtained, since these higher derivatives can be expressed as functions 
of the partial derivatives of/. These would be analogous to higher moments of the spectral 
density function. However, they would be very difficult to evaluate explicitly. 

Apart from explicit evaluations of the expected values of line integrals of this type, the 
rigorous definition of concepts corresponding to spectral density function for the contours, 
regarded as random processes, is a problem of interest. 

4. References (Part 1) 

Longuet-Higgins, M. S. (1952), On the statistical distribution of the heights of sea- waves, J. Marino Res. 9, 

245-266. 
Longuet-Higgins, M. S. (1956), Statistical properties of a moving wave-form, Proc. Cambridge Phil. Soc. 52, 

234-245. 
Longuet-Higgins, M. S. (1957a), Statistical properties of an isotropic random surface, Phil. Trans. Roy. Soc. 

London A250, 157-174. 
Longuet-Higgins, M. S. (1957b), The statistical analysis of a random moving surface, Phil. Trans. Roy. Soc. 

London A249, 321-387. 
Longuet-Higgins, M. S. (1957c), On the velocities of the maxima in a moving wave-form, Proc. Cambridge 

Phil. Soc. 53, 230-233. 
Swerling, P. (1962), Statistical properties of the contours of random surfaces, IRE Trans. Inform. Theory 

IT-8, 315-321. 

Part 2 — Distribution of Specular Points on a Rough 

Spherical Surface 

1. Introduction 

The problem of describing the scattering of high frequency electromagnetic or acoustic 
waves by a statistically rough spherical surface is the primary physical motivation for the 
analysis presented in this article. The theory of scattering from rough surfaces can be traced 
back to the work of Rayleigh [1929], who was interested in the reflection of sound from statis- 
tically rough planes. In comparatively recent times a great deal of attention has been paid 
to the theoretical treatment of similar phenomena involving electromagnetic waves as well 
as acoustical waves, and the subject has been developed to the point where it is now regarded 
as one of the more important aspects of applied scattering theory. 

While many investigators have contributed directly to this subject, both theoretically 
and experimentally, the present article represents an attempt to follow a line of research 
begun by M. S. Longuet-Higgins [1952; 1956; 1957 a, b, c] on a more general class of problems. 
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That work, which has been extended further by P. Swelling [1962], 2 involves the statistical 
description of random surfaces and the resulting geometrical effects. 

Among the questions considered by Longuet-Higgins is the description of those surface 
properties which affect the reflection of coherent unpolarized light according to the theory 
of geometrical optics. He has discussed, for example, the distribution of specular points and 
twinkling effects of a moving sea in relation to certain moments of the surface power spectral 
density [Longuet-Higgins, 1957a]. Although both stationary and moving surfaces were con- 
sidered, the investigations were confined to those cases in which the mean surface is a plane. 

The present article will be concerned with the distribution of specular points for geo- 
metrical optics scattering by a statistically perturbed sphere. The problem will be restricted 
to the case of a source located at a distance which is large compared to the radius of the sphere. 
The procedure followed in the analysis is essentially that of Longuet-Higgins and Swerling. 

An important example of the possible physical significance of the problem treated here 
is in the interpretation of radar scattering data resulting from earth based radar soundings 
of the moon. A complete statistical description of the specular return would also require 
a consideration of the distribution of Gaussian curvature over the lunar surface and, because 
of the coherent nature of the signals used, an analysis of the distribution of the path differ- 
ences associated with the various specular point contributions. However, even the present 
investigation, restricted as it is to the relation between the fundamental properties of the 
perturbed spherical surface geometry and the distribution of specular points alone, should 
provide a certain amount of insight into the nature of the scattering and help to clarify some 
unresolved questions which have been raised concerning the lunar scattering data interpreta- 
tion. In particular, the theory developed here may help eventually to supply the relevant 
connection between the probable lunar surface geometry and the theory of Senior and Siegel 
[1960], in which it has been estimated that between 25 and 30 specular points exist near the 
moon's leading edge. 

2. Surface Angular Correlation Function 

Since the mean surface of the scattering object is a sphere it will be convenient to use 
standard spherical coordinates r, 6, 4> to describe the detailed surface geometry. Thus, if k 
is a fixed unit vector in the direction of the polar angle 0=0, the angle <j> will be measured from 
a fixed radius in the plane through the origin of the coordinate system and perpendicular to k. 
The angle 6 is formed by a radius vector r and a radial line in the direction of k. 

We shall assume that the given closed scattering surface is a random sample from a popu- 
lation of surfaces for which the mean surface is a fixed sphere of radius p, having its center 
at the origin of the coordinate system. 

The population of surfaces can be represented by a stochastic variable of the form 

r(0,<p)=h(P,<p)+ P , (1) 

where h(d, <p) is a random quantity having the mean value zero. 

It will be convenient for obvious geometrical reasons to regard h(6, cp) as represented by 
an expansion in spherical harmonics 

HO, v) = 2 A n P n (cos 6) +S ZJ (Amn sin m<P+B mn cos m<p)P™(cos 6) . (2) 

w=0 n=lm=l 

It will be assumed that the coefficients A mn and B mn are independent random variables having 
zero means: 

E(A mn ) = E(B mn )=0 for all m,n; 

2 Cf. part l. 
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E(A mn A m , n ,) — E(B mn B m , n ,) = for all m,n,m' ,n' except m = m' and^=7i'; 
E(A mn B m , n ,) = for all m,n, m' ,n'; 

in the above relations A 0n and B 0n are to be identified with A n and B n . These assumptions are 
analogous to those which have been made in connection with the planar rough surface scattering 
problem (cf. [Rayleigh, 1929; Longuet-Higgins, 1952, 1956, 1957a; Swerling, 1962; Rice, 1951]), 
wherein a similar formalism based on a Fourier series representation of the perturbations is 
ordinarily used. In the present case the polar singularity of spherical coordinate representations 
dictates the use of the more appropriate spherical harmonics. 

We can define the correlation function relative to h(d, <p) for every pair of surface 
points, determined by the angular coordinates 6, <p, and 0' , <p', as 

oo oo n 

E(h(6, <f,)h(6: <p')) = J2E(Al)P n (cosd)P n (cose')+J2 ^(E(A* mn ) sin mpsin mj 

n=0 7i = l m = l 

+ E(B 2 mn ) cos m^cos m^)^»(cos 0)P£(cos 6'). (3) 

With one additional assumption it will follow that the distribution of h(6 y ip) is isotropic; 
that is, the correlation function E(h(d, <p)h(d' , <p')) is a function only of the angle y between the 
points (0, <p) and (0' </). The angle y satisfies the relation 

cos 7=cos 6 cos 6' +sin 6 sin 6' cos (<p—<p') . (4) 

Because of the uniqueness of the spherical harmonic expansion, the required additional assump- 
tion is 

(n— my 
(n+m) 



E(A* mn ) =E{Bl n ) =2 Z >!;{] W n ; (5) 



the W n defined by (5) are coefficients in the expansion 

E(K(B,<p)W, </>'))=!] W n P n (cosy)=K(cosy)=K(<u), ( M =cos 7 ). (6) 



71 = 



The fact that the relations (5) must hold when the distribution of h(d, <p) is isotropic follows from 
an application of the addition theorem for spherical harmonics to (6) and a comparison of the 
coefficients with those of like terms in (3). 3 

The coefficients W n may be referred to collectively as the spherical spectral power density 
by analogy with the planar rough surface problem. In particular we have from (6) and the 
definition of y an analogous relation for the mean square surface devistion, which we shall refer 
to as the zeroth moment M : 

£W,ri)=SW,=M,. (7) 

n=0 

3. Distribution of Specular Points 

As usual we shall assume that h(6, <p) is a gaussian distributed random variable. The 
probability distribution of h(d, <p) is then determined completely by the quantity M defined in 
(7). Moreover, the joint probability distributions of h(d, <p) and its derivatives with respect to 
6 and <p of various orders are determined in terms of moments of the power spectral density 
function K(fx) defined by (6) (cf. [Swerling, 1962]). Probability distributions for geometrical 
properties of the surface can thus be calculated by the methods described in Swerling's paper. 4 



s The authors are grateful to the referee for pointing out that (6) also implies our other assumptions concerning the correlations of the A m 
and Bmn. 

4 These methods are outlined in part I of this paper. 
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The particular example of the procedure to be considered here is the derivation of the dis- 
tribution of specular points for the scattering of waves by the surface when the source is located 
at a large distance compared to the radius p of the mean sphere. In this connection it will also 
be assumed that p is large compared to the wavelength, so that the dominant scattering mecha- 
nism is geometrical optics, and also that p 2 is large compared to the mean square first and second 
derivatives of h(6, p) with respect to 6. This latter requirement implies that the surface is quasi- 
smooth in some sense. 

We assume that the distant source is located in the direction of the unit polar vector k- 
Then since the vectors Te and r^ are tangent to the surface the quantities 

should be practically zero at a specular point. We then have 

%=he cos 0— (p+A) sin 6, 

rj=h <p cos 0, (8) 

and at a specular point 

£=i?=0. (9) 

We wish to calculate the probability distribution of the event (9). 
First we observe that 

Then following Longuet-Higgins, we consider the joint probability density of the quantities 
£, rj, %e, %<p, r)<p, which we designate by 

This leads at once to the result 

P(Z=v=0)=ddd<p ffJP(0, 0, fo, J„, to)\hi,-$\dhdtjk, ¥ (10) 

for the probability of the occurrence of a specular point in an angular region of measure dddp. 
Now since P(£, 77, £#, £ v „, y^) is a multivariate gaussian distribution it will depend on the 
correlation matrix associated with the five variables, £, rj, &, ^, 77^. The elements of this matrix 
are cross correlations of the variables with 0=0'and <p =<p'. Some of the elements in the 
matrix vanish. 

The correlations can be calculated (cf. [Swerling, 1962]) from (8) and (6). For example, 
we have 

£(£iy)=lim{ E(h e (0, ^M*', *>')) cos cos e' — E(h{B, <p)K.(B' , <p')) sin cos 0' } 
=\im{[K"(n)n 9 Hv'+K' (»)»„>] cos 6 cos e'—K'itinv sin cos 0' } 

=0 
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The matrix, reordered, can be written 

fE(?), E(&,), Efa,), o, 

E(tte), EQS), E{fa 9 ), 0, 

E(£n,), E(6n,), E(n$, 

0, 0, 0, 

v 0, 0, 0, 






0, 

EW), E( V Q 
Eint.) E(£) J 



It follows that P(£, r;, &, £ p , i/ v ) is a product of two simpler distributions: 

Eacli of these is a multivariate gaussian distribution, itself. Then (10) takes the form 



(11) 



; j ({=,=o)=m^ 



;:./:/.■ 



exp {-(a$+2bbrir + c v l+dSl) }\b^^l\db ( K^ (12) 



where the factor A and the exponential coefficients a, b, c, d, are functions of 8 and <p which 
can be obtained from the correlation matrix. 

To evaluate the integral in (12) we must break the region of integration up into regions 
in which (l^^ - O nas a fixed sign. We have then, 

/» CO /* OO /» 00 

exp {— (ax 2 +2bxy-\-cy 2 -\-dz 2 ) }\xy—z 2 \dzdxdy 

«/ — OO tJ — CO «/ — OO 

/» OO /-»CO /» CO 

=4 (z 2 —xy) exp { — ((7z 2 +2&r?/+q/ 2 +<i3 2 ) }dzdxdy 

Jo J -co .7 

+8 f P f^ (Z7/-2 2 ) exp {-(ax 2 +2bxy+cy 2 +dz 2 ) } dzdxdy. (13) 
J o Jo Jo 

The first integral on the right in (13) offers no problem and can be evaluated explicitly. 
We have, in fact, 



and 



Jo J-«Jo 



z 2 exp { - (ax 2 +2bxy+cy 2 +dz 2 ) }dzdxdy= (<jr/d)t(ac-b 2 )-i/8 9 



(14) 



/* CO /* CO /* CO 

xy exp { - (ax 2 +26:n/-f cy 2 +dz 2 ) } dzdxdy = - (6/8) [>ir/(ac-b 2 ) }«-*. (15) 

,/0 J -co JO 

The second integral on the right of (13) can be obtained by differentiating the integral 

/(a, 6, c, r/)= f°° P f Vj " exp {- (ax 2 +2bxy+cy 2 +dz 2 ) } dzdxdy (16) 

Jo Jo Jo 

with respect to b and d. The integral on the right in (16) can be put into a formally simpler 
form. For this purpose, it is convenient to define 

F(s)=dr* erf (dh) = T exp (-dz 2 )dz. (17) 

We now change variables in (16), defining 



u=2xy, 
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(18) 



and obtain 



^=7 I \ F(^uj2) v~ l exp {-[av+bu+cu 2 /(4v)]}dvdu 

4 Jo JO 

=o ^(V^/2) K (u-y/ac) exp (—bu)du. 
* Jo 

In (19) the function .K'o is a Macdonald function, which results from the identity 

r- 1 exp (-aT—Plr)dT=2K (2-yfcJ3), a>0, /3>0. 



(19) 



We are interested in 



~~db 



/i=^t and Tj- 



57 



If we collect the preceding results from (12) through (19) we obtain 

P(Z= v =0)=Addd<p{ ^(r/d)Kac-b 2 )- i -b[T/(ac-b 2 )]id^/2+8r 2 -4T 1 ' 



(20) 



The parameters, a, 6, c, d are functions of 0. They can be expressed in terms of the 
correlation matrix elements. 
Let 

E(?) } E(&e), Efa 9 ) 



and 



Then 



Let 



Ai= 



Em, E(&>, m»i r ) 

E(frt) , E(£»v<e) , E(r)l) 



A; 


I 


EM, 
E(vQ, 


E(vU 
E(itl) 


A=(2tt)-^A^ 




X \E{?), 

fl= 9 


Ar 

E(vl) 1 


i\E(e), 


#(&) 


2 \EfaJ, 


EfatJ 


1 
C= 2 


Em, 


Em 

E(® 


a: 1 


d= l -E{r)A;\ 




M =± 

71 = 


w r 









a: ! 



Mr=T. \n{n-V\)l2}W n 



(21) 



(22) 



M 2 =S {(n-l)n(n+l)(n+2)/8}W n . 

J!=0 
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(23) 



Then we have 

E(F)=M t cos 2 0+(p 2 +M o ) sin 2 0, 

E{&e)=h{p 2 +M,-M x ) sin 20, 

#(&7„) = -(Mi/2) sin 2(9 cos 20, 

S(i 7 2 ) = (M 1 /4)sin 2 20, 

E(v£ v ) = (M 1 /2) sin 2(9 cos 2(9, 

^(^(SMs-Mi+Mo+p 2 ) cos 2 0+4M x , 

S(&^) = (M a /4+M 1 ) sin 2 20, 

S(g) =Afi+ (MV4-M0 sin 2 20, 

#(t7 2 ) = J(M 2 +Af,) sin 2 20. (24) 

We assume that p 2 is large compared to M 0j M u M 2 and that is not near tt/2. Then 

A^M lP 2 sin 2 20 cos 2 0{M 2 cos 2 0+Mj sin 2 0— (M 2 +2M0 sin 2 cos 20}/4, 

A 2 =(M 1 M 2 sin 4 20)/16, 

a~{M x cos 2 0+p sin 2 0) (M 1 +M 2 )-M 2 1 cos 2 20}Z>- 1 (2Af 1 p 2 cos 2 0)" 1 , 

6- - { (Mi cos 2 0+p 2 sin 2 0) (M 2 +4M0 + p 2 M x cos 20}£>- 1 (2M 1 p 2 cos 2 0) "\ 

c-l/(2Z>sin 2 0), (25) 

where 

D=M 2 cos 2 0+M x sin 2 0- (M 2 +2M0 sin 2 cos 20, 

d=2/(M 2 sin 2 20). 

From these results and (21) we find for the limiting case 

p->oo ; 0->O, p sin 0=0(1) 
that 

ac=0(l), 6=0(1), d-^ oo. 
Then 

/f=0(rf-»),/ 1 =0(rf-») J 
and we have 

P( {=1? = 0) ^rf0^{6[7r/(ac-6 2 )]V(2^)+4/}, (26) 

since the other terms are of higher order. Moreover, from (16) we have 



Ii-WW^k 



where 



H= 



exp { — (ax 2 -{-2bxy-{-cy 2 )}dxdy 
Jo Jo 

= (<w?)-* I " j exp {-[x 2 +26(ac)-%+?/ 2 ]}^%. (27) 

Let f==(ac)*, «=6(ac)~*. 

Then 

H=£~ 1 exp { — {x 2 -\-2axy+y 2 )}dxdy 

Jo Jo 

= r 1 (l-a 2 )" i tan- 1 [ (1 -a) /(l +<*)]*, (28) 
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where the integral can be evaluated, for example, by the use of polar coordinates. 
We now have 



I x ~ (1/2) {irld)^{[bl{ac-b 2 )^] tan" 1 [(Vac- ft) /Wac+b)]*- (1/2) (ac-b 2 )' 1 } , (29) 

and therefore 

P(£ =17 =0) (Wir)- 2 Ar H A2' H (2/d)^(ac-b 2 ) -* { irb 

+46 tan" 1 [(^-b)/(^+b)P~2(ac-b 2 y^}ddd<p. (30) 
We have for 0^0: 

A!-MiM 2 P 2 sin 2 0, 

A 2 -MiM 2 sin 4 d, 

a^{{p 2 sin 2 e)(M 2 +M 1 )+M 1 M 2 }l(2M 1 M 2 P 2 ), 

b (2M 2 )- X 

C~(2M 2 sin 2 0)-\ 

d~(2M 2 sm 2 6)- 1 

ac^ { p 2 sin 2 d(M 2 +M 1 ) +MiM 2 }/(±M 1 M 2 2 p 2 sin 2 6) , 

ac-b 2 ~ (p 2 sin 2 6>+M 1 )/(4M 1 Af 2 p 2 sin 2 0). (31) 

It may be observed from (31) and (30) that as 6 approaches zero P(£=7j=0) becomes 
infinite. However, if we ask for the distribution P area of specular points per unit area instead, 
we shall have 

P*re.=Pti=V=0)/(p 2 sin edOdv) (32) 

which approaches a finite constant when 6 approaches zero. 

4. Conclusions 

The distribution of specular points per unit area in the polar region, where 6 is nearly 
zero (although the surface distance pd may be of finite order), is given through expressions 
(30), (31), and (32) in terms of the surface spherical power spectral density moments Mi and M 2 . 
It may be observed by inspection that this distribution becomes indefinitely large if the moment 
M 2 grows indefinitely, and it approaches zero if the moment Mi grows indefinitely. Since M 2 
is related to surface curvature while Mi is related to surface slope these observations are in 
accordance with geometrical intuition. 

The case 0~~ can be treated in a manner similar to the present one, but it is of less interest 

since the distribution will necessarily approach zero rapidly at those angles. 



The authors thank R. Heimiller for his help in several discussions and some concrete 
suggestions which will have contributed materially to any success that this work may enjoy. 
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